Derealization of states in two component superlattices with 

correlated disorder 

T. Hakobyan a , D. Sedrakyan 6 , and A. Sedrakyan a 't I. G6mez c and F. Dommguez-Adame c 
a Yerevan Physics Institute, Br. Alikhanian str. 2, Armenia 
b Yerevan State University 
C GISC, Departamento de Fisica de Materiales, Universidad Complutense, E-28040 Madrid, Spain 

(February 7, 2008) 

Abstract 

Electron and phonon states in two different models of intentionally dis- 
ordered superlattices are studied analytically as well as numerically. The 
localization length is calculated exactly and we found that it diverges for par- 
ticular energies or frequencies, suggesting the existence of delocalized states 
for both electrons and phonons. Numerical calculations for the transmission 

coefficient support the existence of these delocalized states. 
PACS number(s): 73.20.Dx 63.50,+x 71.23.-k 
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I. INTRODUCTION 



Since a remarkable article by Anderson [|TJ, the problem of localization of particles in 
systems with random distribution of parameters is still of continuous interest for physicists. 
It was conjectured by Mott and Twose 0], rigorously proved for some systems [|J and then 
generally argued in Ref. J3] that, in a case of full randomness of the parameters of the model, 
all states are localized in one and two dimensions. However, there exist several exceptions 
to this rule. These exceptions are mainly related to the existence of correlations, either in 
disorder or between the quasi-particles of the system, as well as anomalous (nonexponential) 
localization found at specific regions of the energy spectrum. Recently the interest in the 
investigations of the conditions for breaking of Anderson localization due to correlations 
in the disorder has increased substantially. Evidences were found, that in a presence of 
internal correlations in disordered systems delocalized (extended) states may appear |5|-|2T 



Due to the lack of experimental confirmations, there are some controversies around the 
importance of these results and their physical applications. That is one of the reasons 
why the experimental evidence of extended states, found in the studies of the electronic 
properties of GaAs-AlGaAs superlattice (SL) with intentional correlated disorder by means 



of photoluminescence and vertical d c resistance ||22|| , looks promising. 

Following this line of work, here we consider two component SL's with particular types 
of correlated disorder for thickness of the layers. We demonstrate the appearance of delocal- 
ized states for phonon as well as for electron transport problems. Following the technique 
developed in Ref. |23| we find exactly the transfer matrix for scatterers on the boundaries 
of the layers, and calculate the localization length and dimensionless Landauer resistance, 
which allows us to determine the energy (or frequency) of the resonant states for which 
derealization occurs. Two type of disorder for the thickness of the layers will be considered 
in the paper, namely 

i) The thickness of the one of the SL components (referred to as A layers) is fixed and 
equal to d\, while the thickness of the other component (B layers) is randomly distributed 



2 



with probability 



g{y) = < 



J, 0, < y < di 



(1) 



0, otherwise. 



ii) Again, the thickness of the A layers is set to di, while for B layers we take a sequence 
of fixed and randomly distributed thicknesses. In other words, we take following distribution 
of layers A (fixed) B (fixed) A (fixed) B (random). . . 

In both cases, conditions on energies (and frequencies) of delocalized states are found 
and it is easy to see that they can be fulfilled. We think that these two types of disorder are 
easy to organize in samples grown by molecular beam epitaxy (MBE) and experimentally 
check the existence of extended states for both electrons and phonons, within the spirit of 



Let us consider a SL consisting of two component materials (A and B), grown in the 
x direction with the thicknesses of the layers Axi = Xi — with % = 1, 2, ... , 2N and 
Xj's being the coordinate of the boundaries between the layers. We will investigate the 
propagation of particles and their localization along the grow direction x. The wave equation 
for transversal phonon displacement u(x,y,t) is 



where the velocity of sound c t is defined by the density of matter p and modulus of rigidity 



Solutions of the wave equation (Q) with frequency us are a superposition of forward and 
backward-scattering waves and can be represented as follows: 



Ref. [22 . 



II. TRANSFER MATRIX AND LANDAUER RESISTANCE 



dt 2 



c 2 Au(x,y,t) = 0, 



(2) 



H as c 2 = nj p. 




(3) 




n = 1, . . . , A, 
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where for the wave vector fa , i = 1 , 2 we have 



k 2 -- 



(4) 



Here cu is the velocity of sound in media i; quantities with i = 1 correspond to those of 
material A and, in the same way, quantities with i = 2 correspond to those of material B.In 
equations (HD In (respectively In — 1) numerates the layers B (respectively A). 

We should now impose the boundary conditions on the solutions (|3|) for displacements 



U2n and U2 n -i, as it was demonstrated in Ref. |23 



n—l, 



flld x U2n-l(x2n-l) , 



U2n{x2n-l) ~ U2n-1 {%2n-l / 



n=l,2,---N, 



(5) 



which are nothing but continuity conditions on the displacements Ui(x,y,t) and the forces 
Hid x Ui(x,y,t) at the boundaries of the layers. 

The solution of the boundary conditions allows us to express linearly all amplitudes 
Cj and Q of scattering modes in the slice i through amplitudes of the initial wave c\ and C\ 
as follows 



^2n = II T i^l = T ^l' 

i=i 



(6) 



where we have defined 



i = 1, 2, . . . , 2N. 



(7) 



The expression of the transfer matrix Tj 

Tj = T 2j T2j-i 



(8) 



with 



a 



-( + —j^) Sm fcl(^2i-l - %2j-2) + COS k l {x 2j -l - X2J-2) 



D ik2(x 2 j-2-X2j-3) 



(9) 



is easy to obtain performing the product of T 2 j and T 2 j-\ matrices 



T, 



1 



2.7 



T, 



2fJ, 2 k 2 



2/jiki 



(juifci + /i 2 fc 2 )e ifcl( ^- 1 -^- 2) (/i 2 fc 2 - /iiJfeOe-** 1 ^- 1 - 3 **- 3 ) 



v (M2 - /ii^Oe^ 1 ^- 1 -^- 2 ) (//life! + /i2 A; 2 ) e - ifc i(^i-i-^-2) 



(//ifci + fA 2 k 2 )e ik ^ X2 i- 2 - x v- 3) (frfa - fi 2 k 2 )e~ ik ^- 2 - X2 i-^ 



(10) 



If we now focus on a model where electrons with effective mass m ; and potential energy 
Vi at layer A(i = 1) or B(J = 2) impinge on the SL, then it is necessary to change the 
equation of motion for phonons by the Schrodinger equation for the envelope function 
(h = 1 and k± = hereafter): 

,du(x,t) f 1 



dt 



K . lm &-V i )u{x,t) = U 1 , 



1,2,...,2A, 



(11) 



but the form of general solution (RJ) is valid provided 



k] =2m i (E-V i ). 



Analogously, the boundary conditions (||) now will read as 



(12) 



-d x u 2n _i(x 2n _i), n = l,2, 



N. 



(13) 



u 2n {x 2n -i) — U 2n -.\\X 2n -i) 
—dxU2n{x2n-l) = — 

m 2 rrii 

It is a matter of simple algebra to see from equations fll2|) and ([U|) that the transfer ma- 
trix (P) for the electron problem will have the same form as for phonons in the expression (|8| 
but replacing fii — > l/rrii. 



In Ref. |2B] the problem of transport of particles in the one dimensional space for a wide 
class of disordered models was considered within the transfer matrix approach and general 
results were obtained. It was demonstrated that the transfer matrix of one dimensional 
problems belongs to SL(2, R) group and randomness of media can be exactly taken into 
account for such quantities as the Landauer resistance f24|. It is easy to see that our 
transfer matrix also belongs to SL(2, R) group. 



Following Refs. [ 24|— 126 1 let us now define the dimensionless Landauer resistance as the 
ratio of reflection and transmission coefficients 



P 



T 



12^21- 



(14) 



In order to proceed further we should use the relation for the direct product T ® T 1 
presented in Ref. [^3j 



(15) 



where 



(16) 



is the spin-one part of the direct product. But for Landauer resistance we need to calculate 
T g) T + . It is easy to see from (||) that 



Therefore, by multiplying flT5|) from the left and right by cr 3 we will have 

1 



(17) 



(18) 



It is straightforward now to calculate the Landauer resistance p by using formulas (|T^)-(|T8|). 
which seem to depend only on (3, 3) element of the product of the transfer matrices 



P= 2 



N 



-l + CQA, 



(19) 



This expression is of remarkable interest because it is multiplicative in Aj U (j = 1,2,..., N), 
each of which depends only on local parameters (thickness and the other model parameters) 
of the j th pair of layers. Therefore, this expression for the Landauer resistance becomes 
valid for media with arbitrary distribution of the parameters. 
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III. CONDITIONS FOR THE EXISTENCE OF EXTENDED STATES 



As we see from the expression (|19"D , due to the multiplicative form of the dependence 
of the Landauer resistance on local parameters, its average over different type of correlated 
disorder can be easily taken. We should simply average Af u in each layer separately and 
then take the product of them. Now we will consider the two type of correlated disorder 
mentioned in the Introduction. 

Let us take fixed thickness for the A layers of the SL as d\ and a random distribu- 
tion of thicknesses for the component B by use of the probability distribution (Q). Simple 
substitution of expression (||) for the transfer matrix T into the formula ([16]) for Aj U gives 

A) 1 = cos[2/M{] cos[2Jfe 2 d£] - cosh[0] sm[2k 2 4] smpJ^dj], 
Af = (cosh 2 [#] cos[2M{] - sinh 2 [fl]) sin[2fc 2 4] + cosh[#] cos[2M 2 ] sin[2Mi], 
A] 3 = -zsinh[#] [cos[2Jfe 2 4] sin^Mi] - cosh[0](l - cos[2M{]) sin[2fc 2 4]] , 
Af = - sin[2fc 2 4] cos[2Mi] - cosh[#] cos[2A; 2 4] sin[2Mi], 
Af = (cosh 2 [fl] cos[2M{] - sinh 2 [fl]) cos[2Jfe 2 4] - cosh[#] sin[2A;i4] sin[2A; 2 4], 
Af = ismh{6] [sin^Mi] sin[2fc 2 4] + cosh[0](l - cospMl]) 
Af = ismh[6] sin^/Ml], 
Af = isinh[20](l - cos[2k 1 d{])/2, 

Af = cosh 2 [0] - sinh 2 [0] cos[2k x d{l (20) 
where 6 is defined by 

cosh[0] = 1 ( ^ + ^l). (21) 

and d\ = x 2 j-i — %2j-2 and d 2 = x 2j _ 2 — x 2 j_ 3 are the thicknesses of the j th pair of layers. 

Now we should fix d\ = d\ for the component A and take the average over d 2 using the 
probability distribution ([I]), which will give (A)^ defined by the expressions (^), but where 
cos[2k 2 d 2 ] and sin[2/c 2 <i 2 ] are changed by their average values 



/ ro; Jin sm[2d 2 k 2 } 



2g? 2 A; 2 

6 = (sin[2fc 2 d 3 2 ]> = , , - ■ (22) 
Then, for the averaged Landauer resistance we will have 

( P ) = \[- i + m N r] ■ (23) 

For large sample sizes (N ^> 1), as it was first argued in Refs. |2"3] , |2"7|] , the resistance 
should behave as e lN , where the Lyapunov exponent 7 provides the phonon correlation 



length. Using fl23[) and the definition of Lyapunov exponent 7 = liniAr^oo In p/N we can find 



an exact expression for localization length 

T^lnA, (24) 

where A is the closest to one eigenvalue of the matrix (A 33 ). Excitations are localized or 
not depending on the behavior of £. If at some frequency uj c (or energy) the localization 
length becomes infinite, we generally have delocalized states |28| and the expression 



shows that it will occur when X(u c ) = 1. Therefore we should elucidate whether the matrix 
(A) can support unity eigenvalue or not. It is then necessary to calculate the determinant 
of the matrix 1 — (A) 



det [1 - (A)] = - sin 2 ^) ( ^ - (a 2 + 6 2 - 1 

2 \p2k2 HihJ v 



(25) 



from where it follows that the condition to have an extended state is 

sin(d x A;i) = 0. (26) 

Let us now fix the thickness of the component A of the SL as d\ and for the component 
B take fixed and random thicknesses in a sequence. Then the extended states can appear 
when 

det[l - A(A)] = 0, (27) 
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where A is matrix ( pO|) with fixed thickness of B layers d 2 = d 2 . It turns out that the 
condition ([27D is equivalent to the equation 

cos kidi cos k 2 d 2 ( ^ 1 1 + ^ 2 2 | sin fci<ii sin fc 2 ^2 = 0. (28) 

IV. NUMERICAL RESULTS 

In order to validate the results of our previous formalism, we performed some numerical 
calculations which allow us to show the existence of the extended states discussed above. 
We will focus our attention on the electronic model of disorder we previously referred to as 
model i), and calculate for that kind of disorder the transmission coefficient as a function 
of energy, as well as a function of the system size when the energy is fixed to one of that 
given by expression fl26|). The transmission coefficient was numerically computed using the 



transfer matrix formalism 29,30 



Fig. [I] shows the transmission coefficient calculated in model i) as a function of the 
reduced energy E/V 2 for states above the barrier. We have chosen a GaAs-AlGaAs SL as a 
typical example with the following structural parameters: d\ = 200A, d 2 = 15A, V 2 = 0.4eV 



and A" = 200. The arrows are located at the energies predicted by relation (|26|). It turns 
out that these energies are given by 

n 2 7r 2 h 2 , , 

E =2^T (29) 

n being an integer number. It is clear that they coincide with the sharp resonances in the 
transmission coefficient that can be observed in the figure. 

To check whether the energies given by the previous relation (^9]) correspond to extended 
states or not, we represent in Fig. [2] the transmission coefficient for a couple of such energies 
as a function of the size of the system, and compare it with the case in which the energy 
of the state lies between two of them. For the energies in (|29| ) the transmission coefficient 
remains constant as a function of the size N, this behavior being expected for an extended 
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state. Meanwhile, for a state with energy between two resonances the transmission coefficient 
decays exponentially. 



V. CONCLUSIONS 

In summary, we have shown that two particular models of correlated disordered SL's 
exhibit delocalized states for electrons as well as for phonons. This result has been demon- 
strated analytically as well as numerically. We have found exactly the energy and frequency 
for which extended states appear. Notice from Fig. [I] that the resonances of the transmission 
coefficient around the theoretical values fl26|) are rather broad. This suggests that electron 
and phonon states close to the values given by (|2~6"D should display a rather large localization 
length, even larger than the SL length. This is relevant for transport measurements provide 
the Fermi level (in the case of electrons) are located close to one of these maxima. In such 
a case, one would expect an enhancement of the dc conductance of the sample, as it was 



actually observed in the case of the so-called random dimer SL's [22 
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FIGURES 

FIG. 1. Transmission coefficient, r, as a function of the reduced energy EjVi- The arrows 



show the energies given in (29). The height and the nominal width of the barriers are, respectively, 
Vi = 0.4 eV and di = 15 A, the width of the wells is d\ = 200 A, and the number of periods is 
N = 200. 

FIG. 2. Transmission coefficient, r, as a function of the system size N . The dotted and dashed 
lines correspond, respectively, to n = 7 and n = 14 in (|29|). The solid line corresponds to an energy 
between n = 14 and n = 15. Structural parameters are the same than in the previous figure. 
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